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The correlation length of the square- lattice spin-1/2 Heisenberg antiferromagnet is studied in the low- 
temperature (asymptotic-scaling) regime. Our novel approach combines a very efficient loop cluster 
algorithm - operating directly in the Euclidean time continuum - with finite-size scaling. This 
enables us to probe correlation lengths up to ^ « 350, 000 lattice spacings - more than three orders 
of magnitude larger than any previous study. We resolve a conundrum concerning the applicability 
of asymptotic-scaling formulae to experimentally- and numerically-determined correlation lengths, 
and arrive at a very precise determination of the low-energy observables. Our results have direct 
implications for the zero-temperature behavior of spin-1/2 ladders. 
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Soon after the discovery of high-temperature supercon- 
ductivity in doped lamellar copper oxides it was found 
that the undoped compounds are quasi-two-dimensional 
(2-d) spin S* = 1/2 quantum antiferromagnets. A theo- 
retical model that captures the essential features of these 
materials is the nearest-neighbor quantum antiferromag- 
netic Heisenberg model (AFHM) on a square lattice. 
Through experimental, numerical, and theoretical efforts 
much progress has been made in the understanding of 
these systems. In particular, detailed neutron scattering 
measurements of the spin-spin correlation length in the 
magnet Sr2Cu02Cl2 were found to be described quantita- 
tively ||l| by both high-temperature numerical results for 
the AFHM |Q and low-temperature theory for the renor- 
malized classical regime of the (2 + l)-d 0(3)-symmetric 
non-linear cr- model [^|| . 

The ground state of the above systems shows long- 
range antiferromagnetic order, thus spontaneously break- 
ing the 0(3) rotational symmetry to 0(2). The 
low-energy excitations are two massless bosons called 
magnons or spin-waves. These long-range excitations 
determine the dynamics at low energies. One can use 
chiral perturbation theory (CPT) to derive universal ex- 
pressions for low-energy observables in terms of three 
material-specific parameters: the staggered magnetiza- 
tion A4s, the spin- wave velocity c, and the spin stiffness 
Ps Both numerical data and predictions of CPT are 
in apparent good agreement with experimental results 
for S = 1/2 01 . However, neutron scattering measure- 
ments on antiferromagnets with 5* > 1/2 [^]j6j reveal a 
striking discrepancy with CPT predictions based on 3- 
loop asymptotic scaling. It has been suggested that for 
S > 1/2, asymptotic scaling sets in only at very low tem- 
peratures - that is, for correlation lengths much larger 
than those accessed experimentally and numerically 0. 

In this Letter, we investigate the correlation length of 
the 2-d nearest-neighbor square-lattice spin-1/2 Heisen- 
berg antiferromagnet at unprecedentedly low tempera- 



tures. We resolve the puzzle concerning the applicability 
of asymptotic scaling to experimentally- and numerically- 
determined correlation lengths [^||^-|^. To effect this, we 
combine an efficient loop cluster algorithm operating in 
the Euclidean time continuum - hence with zero system- 
atic error - with a finite-size scaling technique. With 
this powerful new approach, the infinite-volume correla- 
tion length ^ is probed up to « 350, 000 lattice spacings 
a. This allows us to paint a detailed picture of the varia- 
tion of ^ with temperature well into the asymptotic scal- 
ing regime, and to arrive at the presently most precise 
determination of the low-energy observables of the spin- 
1/2 AFHM. We find that asymptotic scaling of with 
the 3-loop /3-function of the 2-d classical 0(3) model sets 
in only at about lO^a, while a 4-loop fit describes the 
data already at ^/a w 100 — 200, corresponding to the 
largest correlation lengths measured experimentally [0. 
By exchanging a spatial with the Euclidean time direc- 
tion, as suggested in Ref. |9|, our results may be appHed 
to spin-1/2 ladders at zero temperature. 

In CPT, the 2-d quantum spin system is described 
by a (2 -I- l)-d 0(3) symmetric Euclidean field theory. 
At non-zero temperature T the Euclidean time direction 
has a finite extent 1/T. For a system of massless par- 
ticles - that is, one with an infinite zero-temperature 
correlation length - the non-zero temperature system 
appears dimensionally reduced to two dimensions, be- 
cause 1/T is then negligible compared to ^. However, 
the Hohenberg-Mermin-Wagner-Coleman theorem for- 
bids interaction among massless Goldstone bosons in two 
dimensions. Consistent with this, the 2-d 0(3) model is 
known to have a non-perturbatively generated mass gap. 
This in turn implies that the spin- waves of a 2-d quantum 
antiferromagnet at non-zero temperature also acquire a 
mass. Using perturbative renormalization group argu- 
ments, Chakravarty, Halperin, and Nelson (CHN) |^ de- 
rived the expression 
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Consistent with dimensional reduction, ^ is exponentially 
large compared to c/T. We note that the exponent in 
Eq.(|l]) comes from a 1-loop calculation, while the factor 
cjlups is a 2-loop result. 

Hasenfratz and Niedermayer Q averaged the (2 + 1)- 
d field over cubic space-time volumes of size l/T in the 
Euclidean time direction and c/T in the two spatial di- 
rections. Since at low temperatures ^ ^ c/T, the field 
is essentially constant over these blocks. The averaged 
field naturally lives at the block centers, which form a 2- 
d lattice of spacing a! — c/T (which is different from the 
lattice spacing a). Hence, the effective action of the aver- 
aged field defines a 2-d lattice 0(3) model. Using CPT, 
as well as the exact mass-gap and the 3-loop /3-function 
of the 2-d 0(3) model Hasenfratz and Niedermayer 
extended the CHN-formula to 
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which we call the CH2N2-formula. This equation is 
valid at low temperatures, that is, for large correlation 
lengths. When the correlation length is correctly de- 
scribed by Eq.(||), it scales asymptotically with the 3-loop 
/3-function of the 2-d 0(3) model. The undetermined 
OiT"^ / p1) term represents a 4- loop effect. 

In relativistic quantum field theory the lattice spacing 
serves as an ultraviolet cut-off that is ultimately removed. 
The question of asymptotic scaling is thus unphysical be- 
cause it involves the bare coupling constant. In the 2-d 
AFHM, however, the lattice spacing of the induced clas- 
sical model is a' = c/T ^ where T is the physical tempera- 
ture. Hence, the question of asymptotic scaling becomes 
physical. A priori, it is unclear for what values of f one 
should expect asymptotic scaling for the effective 2-d lat- 
tice action. 

This raises the important question: for what ^ does the 
CH2N2-formula work correctly? We address this issue for 
the AFHM, defined by the Hamiltonian 
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(3) 



where J > is the antiferromagnetic coupling, Sx is a 
quantum spin- 1/2 operator located at point a; of a square 
lattice with spacing a, and jl is the unit vector in the /x- 
direction. Previously, comparisons of Monte Carlo calcu- 
lations of the staggered and uniform susceptibilities {xs 
and Xu) with CPT predictions yielded c = 1.68(l)Ja, 
Ps = 0.186(4) J, andA^, = 0.3083(2)/a2 [|ll|,0. The cal- 
culation of Rcf. [0 indicates that CPT results for Xs and 
Xu of the same kind as Eq.(||) - but with known 0{T'^ / pi) 
terms - work very weU for the AFHM for T < 0.2J. 
Hence, one expects that asymptotic scaling at the 4-loop 



level sets in at about the same temperature. At T = 0.2 J 
the CH2N2-formula predicts ^/a « 150. To avoid finite- 
size effects in numerical simulations, one must simulate 
lattices with a spatial extent at least L « 6^. Even 
with the best algorithms available |l^,|ll|,^, working 
on lattices starting with = 900^ is extremely time- 
consuming. 

How then can we investigate ^ in the low-temperature 
regime? We make use of finite-size effects instead of try- 
ing to avoid them. Finite-size scaling methods were used 
by Kim |Q and by Caracciolo et al. |^ to test asymp- 
totic scaling of the correlation length with the bare cou- 
pling constant in the classical 2-d lattice 0(3) model. To 
investigate ^ in the AFHM, we apply the technique of 
Ref. |15| , in which renormalization group methods are 
used to justify the search for a universal scaling function 



^(2i)/^(T) = T(e(T)/L). 
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Carracciolo et al. |15| showed that their data indeed col- 
lapse to a single scaling function F{£^{L)/ L). They also 
showed that iteration of Eq. (^ yields rapid convergence 
to ^, the correlation length of the infinite system. 

The universal function F for the classical 2-d 0(3) 
model was determined very precisely in Ref. . Since at 
low temperatures the quantum model reduces to a classi- 
cal 2-d lattice 0(3) model, one can use the same function 
to deduce ^ from £,{L) data. We verified that indeed the 
same scaling function works for the Heisenberg model 
with measurements of £,{L) for lattices of successively 
larger size. This calculation yields very good, though 
not exact, agreement with F{£^{L)/ L). We incorporated 
a tiny correction to the scaling function to account for 
these violations at each stage of the iteration. We found 
that the correlation length so determined was insensitive 
to the form of the scaling violation, for several different 
fitted functional forms. 

We emphasize that Eq.(^ assumes universal behavior, 
that is, scaling, but not asymptotic scaling. This is im- 
portant, because we want to use Eq.(^ to compare with 
the CH2N2-formula Eq.(^ without bias. 

The scaling procedure is very sensitive to small changes 
in the finite- volume correlation length ^(T), that is, a 
small error in ^(T) can lead to large uncertainties in the 
infinite- volume ^. Hence, one needs a very accurate nu- 
merical method to determine £,{L). Fortunately, for the 
AFHM there exists a very efficient loop cluster algorithm 
|]l3| , p| , which practically eliminates auto-correlations in 
successive Monte Carlo configurations. The cluster algo- 
rithm also enables improved estimators which drastically 
reduce statistical errors. Finally, implementing the clus- 
ter algorithm in continuous Euclidean time completely 
eliminates the systematic error due to the Trotter- Suzuki 
discretization of time [T^ . This continuous-time cluster 
algorithm (CTCA) has greatly reduced storage and com- 
puter time requirements, enabling simulation at very low 
temperatures. 
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For J/T = 0.5, 1.0, 3.5, we carry out simulations 
that measure the correlation length directly in a large 
volume L k, 6^, while for J/T = 4.0, 4.5, 12.0, we 
use the finite-size scaling technique. A single-cluster ver- 
sion of the CTCA is used with an improved estimator 
for the staggered correlation function, then £,{L) is ex- 
tracted from the correlation function using the second- 
moment method [ p^|JT5| . For the direct-measurement 
cases, 10^ Monte Carlo configurations are used; for the 
scaled-measurement cases, 4 x 10^ are used. The nu- 
merical data for ^(L) together with the inferred infinite- 
volume ^ are given in Table ^, and are compared with 
experimental data and the CH2N2-formula in Fig.|l|. The 
largest accessible ^ is more than three orders of magni- 
tude larger than those of any previous study 0,§,|l6|,|| . 

Deviations from asymptotic scaling are invisible on the 
semi- log scale of FigJ^. Figure || shows the deviation from 
2-loop asymptotic scaling as a function of T/2'Kps. Here, 
the 3-loop result of Eq.(||) is a line of slope —1/2. The 
quantitative agreement shown in Fig.^ between experi- 
ment and theory in the regime 2.4 < J/T < 5.3 is seen 
to be coincidental. We find that with decreasing tem- 
perature the correlation length crosses and undershoots 
the 3-loop result before approaching it from below. Con- 
sequently, the deviations at intermediate and high tem- 
peratures from the CH2N2-formula turn out to be rel- 
atively small. Although neutron scattering results for 
Sr2Cu02Cl2 cover a wide range of ^ with relatively small 
errors 0|, the deviations from 3-loop asymptotic scaling 
turn out to be too subtle to be discernable experimen- 
tally. This explains why the low-temperature result of 
Eq. appears to describe the experimental data at such 
high temperatures. 

Owing to the exponential dependence on ps , the place- 
ment of the numerical data on the graph of Fig.^ is ex- 
tremely sensitive to ps in the low-T regime. This enables 
a very precise estimate of ps . The correlation length data 
are fitted simultaneously with previously obtained results 
for Xs and Xu in cubical |jll| and cylindrical jl^ space- 
time geometries, resulting in 

c = 1.657(2) Ja, ^ 0.1800(5) J, Ms = 0.30797(3)/a2. 

This is in very good agreement with results of an expan- 
sion around the Ising limit, which gives c = 1.654(11) Ja, 
Ps = 0.182(5)J, and Ms = 0.307(l)/a2 0. 

In order to fit the correlation length data down to 
f/a w 100, it was necessary to include the 0{T^/pl) 
4-loop term - not determined in Ref. - with a fit- 
ted coefficient C2. We find C2 = -0.75(5) if the fit is 
restricted to quadratic order and to J/T > 2.8. Higher- 
order fits (including C3, C4, ...) shift the value deter- 
mined for C2, although it remains 0(1) and stable. A 
previous study also found that quadratic-order terms are 
needed to match the precision of the Xs and Xu mea- 
sured with the CTCA [p2|. Thus, it should come as no 



surprise that asymptotic scaling at the 3-loop level of the 
CH2N2-formula sets in at rather large correlation lengths 
of ^ w lO^a. The same is true for the classical 2-d 0(3) 
model with the standard lattice action |l^ . 

We also carried out an additional independent fit for 
the prefactor in the CH2N2-formula. That is, we multi- 
plied Eq.(^ by an extra fitted parameter A, and found 
A — 1.02(2) (using the 4-loop form for the /3-function) , 
in excellent agreement with Eq.(^). This gives added 
confidence both in our scaling procedure and in the cor- 
rectness of the exact massgap. 

A large discrepancy between the CH2N2-formula and 
experimental data has been discovered for systems with 
5 > 1/2 A likely explanation of this discrepancy is 

that 3-loop asymptotic scaling again sets in only at very 
small temperatures. In fact, the existence of a crossover 
from high- to low-temperature behavior was suggested in 
Ref. . Although the asymptotic region is inaccessible to 
experiments, we have demonstrated that it is possible to 
investigate it using the cluster algorithm in combination 
with finite-size scaling. A similar study for S > 1/2 is 
presently in progress. 

Recently, it was pointed out that the CH2N2-formula 
can also be applied to the correlation length at T = 
of even- width Heisenberg spin-ladders ||]. These ladders 
comprise an even number n of coupled S = 1/2 chains 
with periodic boundary conditions. This mapping ex- 
changes one spatial direction with the Euclidean time 
direction, so the time extent plays the role of n. The 

3- loop /3-function is the same for the square-lattice and 
the ladder. In contrast, the undetermined higher-order 
terms are infiuenced by the distinction that space is a 
lattice, whereas Euclidean time is continuous. Hence, 
the 0{T'^ / p1) terms are expected to be different for the 
square-lattice and the spin-ladders. Early density-matrix 
renormalization group calculations of the spin gap of 
even- width ladders were reported in Ref. [^s) . The corre- 
lation length of spin-1/2 ladders was studied numerically 
for n < 6 [^^|. We see in Fig.|| that 3-loop asymp- 
totic scaling sets in only at T « 0.13 J. We thus expect 
the 3-loop CH2N2-formula to work quantitatively only 
for ladders with widths greater than n = c/T « 12. It 
was already realized in Ref. ||l^ that the 3-loop result is 
insufficient to describe the numerical data for n < 6. A 

4- loop form should work for n > 6. 

In conclusion, we have combined a powerful and ac- 
curate quantum Monte Carlo technique with finite-size 
scaling. This has made possible the determination of 
the correlation length of the S* = 1/2 nearest-neighbor 
square-lattice antifcrromagnctic Heisenberg model at un- 
precedentedly low temperatures. Our study resolves the 
conundrum concerning the applicability of the 3-loop 
asymptotic scaling description for S — 1/2, and it points 
the way toward a quantitative resolution of this issue for 
S > 1/2. Our results also have direct implications for 
the low-temperature properties oi S = 1/2 ladders. 



3 



Wc arc indebted to A. Ferrando and F. Niedermayer 
for very interesting discussions. We also thank A. Pelis- 
setto for providing us with the universal finite-size-scaling 
function F. This work was supported by the NSF un- 
der Grant No. DMR 97-04532, the International Joint 
Research Program of NEDO (New Energy Development 
Organization) International Joint Research Grant, and 
the DOE under cooperative research agreement number 
DE-FC02-94ER40818. U.-J.W. is also supported by an 
A. P. Sloan fellowship. 



1.5 


80 


1.818(4) 


1.818(4) 


2.0 


80 


3.351(8) 


3.351(8) 


2.5 


80 


6.23(1) 


6.23(1) 


3.0 


160 


11.60(3) 


11.60(3) 


3.5 


160 


21.23(2) 


21.25(2) 


4.0 


80 


31.3(1) 


39.4(2) 


4.5 


80 


39.0(1) 


69.9(7) 


5.0 


80 


44.6(1) 


126 (2) 


6.0 


80 


52.5(1) 


403 (6) 


7.0 


80 


58.8(1) 


1275 (27) 


8.0 


80 


64.2(1) 


3960 (120) 


10.0 


80 


73.5(1) 


38460 (1010) 


12.0 


160 


159.0(2) 


352000 (10500) 



[1] M. Greven et al, Phys. Rev. Lett. 72, 1096 (1994); Z. 

Phys. B 96, 465 (1995). 
[2] M.S. Makivic and H.-Q. Ding, Phys. Rev. B 43, 3562 

(1991). 

[3] S. Chakravarty, B.I. Halperin, and D.R. Nelson, Phys. 

Rev. B 39, 2344 (1989). 
[4] P. Hasenfratz and F. Niedermayer, Phys. Lett. B 268, 

231 (1991). 

[5] P. Hasenfratz and H. Lcutwylcr, Nucl. Phys. B 343, 241 
(1990); P. Hasenfratz and F. Niedermayer, Z. Phys. B 
92, 91 (1993). 

[6] K. Nakajima et al, Z. Phys. B 96, 479 (1995); Y.S. Lee 

et al., to be published. 
[7] N. Elstner et al, Phys. Rev. Lett. 75, 938 (1995). 
[8] J.-K. Kim, D.P. Landau, and M. Troyer, Phys. Rev. Lett. 

79, 1583 (1997). 
[9] S. Chakravarty Phys. Rev. Lett. 77, 4446 (1996). 
[10] P. Hasenfratz, M. Maggiore, and F. Niedermayer, Phys. 
Lett. B 245, 522 (1990); P. Hasenfratz and F. Nieder- 
mayer, Phys. Lett. B 245, 529 (1990). 
[11] U.-J. Wiese and H.-P. Ying, Z. Phys. B 93, 147 (1994). 
[12] B.B. Beard and U.-J. Wiese, Phys. Rev. Lett. 77, 5130 
(1996). 

[13] H.G. Evertz, G. Lana, and M. Marcu, Phys. Rev. Lett. 
70, 875 (1993). 

[14] J. K. Kim, Phys. Rev. Lett. 70, 1735 (1993); Phys. Rev. 

D 50, 4663 (1994). 
[15] S. Caracciolo et al, Phys. Rev. Lett. 75, 1891 (1995). 
[16] M. Greven, U.-J. Wiese, and R.J. Birgeneau, Phys. Rev. 

Lett. 77, 1865 (1996). 
[17] Z. Weihong, J. Oitmaa, and C.J. Hamer, Phys. Rev. B 

43, 8321 (1991); C. Hamer, Z. Weihong, and J. Oitmaa, 

Phys. Rev. B 50, 6877 (1994). 
[18] S.R. White, R.M. Noack, and D.J. Scalapino, Phys. Rev. 

Lett. 73, 886 (1994). 
[19] O.F. Syljuasen, S. Chakravarty, and M. Greven, Phys. 

Rev. Lett. 78, 4115 (1997). 

TABLE I. Finite- and infinite-volume correlation length. 
Infinite- volume ^/a are deduced using finite-size scaling. 

J/T L/a i.{L)/a ^/as^(oo)/a 

0.5 20 0.481(1) 0.481(1) 

1.0 80 0.973(2) 0.973(2) 
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FIG. 1. Experimental [g| (filled circles) and computed val- 
ues (open circles; errors are much smaller than symbol size) 
of the staggered- spin correlation length. The solid line is the 

3- loop CH2N2-formula, Eq.(||), with = 0.1800. 

FIG. 2. Analysis of the correlation length as the asymp- 
totic scaling regime is approached. To remove the exponential 
inverse-temperature dependence, ^ is divided by the 2-loop 
formula, using the best-fit value ps — 0.1800(5). The 3-loop 
result (straight line) is accurate only for ^/a « 10^, while the 

4- loop regime (parabola) begins at roughly ^/a « 10^. 
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FIG. 1. Experimental [1] (filled circles) and computed values (open circles; errors are much smaller than 
symbol size) of the staggered- spin correlation length. The dashed line is the 3-loop CH2N2 formula, Eq.(2), 
with p,=0. 1800. 
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FIG. 2. Analysis of the correlation length as the asymptotic scaling regime is approached. To remove the 
exponential inverse-temperature dependence, ^ is divided by the 2-loop formula, using the best-fit value 
yOs=0. 1800(5). The 3-loop result (straight line) is accurate only for ^a«10^, while the 4-loop regime (parabola) 
begins at roughly =10^. 
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